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^ : Abstract 



We consider a five dimensional DGP-brane scenario endowed with a non-minimally coupled 
scalar field within the context of Brans-Dicke theory. This theory predicts that the mass appearing 
in the gravitational potential is modified by the addition of the mass of the effective intrinsic 
' curvature on the brane. We also derive the effective four dimensional field equations on a 3 + 1 

dimensional brane where the fifth dimension is assumed to have an orbifold symmetry. Finally, we 
discuss the cosmological implications of this setup, predicting an accelerated expanding universe 
O , . with a value of the Brans-Dicke parameter ui consistent with values resulting from the solar system 



observations. 



1 Introduction 



[ Over the past decade the possibility of the observable Universe being a brane-world [1] embedded 

I in a higher dimensional space-time has generated a great amount of interest. This was motivated 

OO ■ by the fact that there is a strongly coupled sector of x Eg heterotic string theory which can be 

described by a field theory hving in 11-dimensional space-time [2]. The 11-dimensional world consists 
of two 10-dimensional hypersurfaces embedded on the fixed points of an orbifold and the matter 
fields are assumed to be confined and live on these hypersurfaces which are known to be 9-branes. 
. After compactification of the 11-dimensional theory on a Calabi-Yau 3- fold, we obtain an effective 

5-dimensional theory [3] which has the structure of two 3-branes located on the orbifold boundaries. 
This scenario has motivated intense efforts to understand the case where the bulk is a 5-dimensional 
anti de-Sitter space. In this setup, gravitons are allowed to penetrate into the bulk but are localized 
on and around the brane [4] . It was then shown that in a background of a non-factorizable geometry 
an exponential warp factor emerges which multiplies the Poincare invariant 3-1-1 dimensions in the 
metric. The model consists of two 3-branes situated along the 5th dimension, compactified on a 
/Z2 orbifold symmetry where the two branes must have opposite tensions. The evolution equation 
followed from such a brane scenario differs from that of the standard four dimensional evolution 
equation when no branes are present [6]. The existence of branes and the requirement that matter 
fields should be localized on the brane lead to a non-conventional cosmology, necessitating a more 
concerted study. A large number of studies have been devoted to the effective gravity induced on the 
brane [7] and, in particular, a great amount of interest was generated on inflationary cosmology [8]. 
More recently, post-inflationary brane cosmology has been also considered in [9]. Not surprisingly, 
the problem of the cosmological constant has become a focal point in the brane-world studies where, 
for example, in [10, 11, 12] a flve dimensional action with a scalar fleld is non-minimally coupled 
to flve dimensional gravity and to the four dimensional brane tension. There has also been some 
discussion on the localization of gravity [13]. A feature common to these type of models is that they 
predict deviations from the usual 4D gravity at short distances. 
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A somewhat different approach within the brane- world framework is the model proposed by Dvali, 
Gabadadze and Porrati (DGP) [14, 15]. It predicts deviations from the standard AD gravity over 
large distances. The transition between four and higher-dimensional gravitational potentials in the 
DGP model arises because of the presence of both the brane and bulk Einstein terms in the action. 
The Friedmann-like equations governing the cosmological evolution of a brane possessing an intrinsic 
curvature term in its action have already been derived and discussed for an AdS-Schwarzschild bulk 
space-time [16]. Cosmological consideration of the DGP model was first discussed in [17] where it was 
shown that in a Minkowski bulk space time wc can obtain self-accelerating solutions. In the original 
DGP model it is known that 4Z) general relativity is not recovered at linearized level. However, some 
authors have shown that at short distances we can recover the 41? general relativity in a spherically 
symmetric configuration, see for example [18]. An important observation was made in [19, 20] where 
it was shown that the DGP model allows for an embedding of the standard Friedmann cosmology 
whereby the cosmological evolution of the background metric on the brane can entirely be described by 
the standard Friedmann equation plus energy conservation on the brane. This was later extended to 
arbitrary number of transverse dimensions in [21]. For a comprehensive review of the phenomenology 
of DGP cosmology, the reader is referred to [22] . 

It is worth mentioning that an interesting feature of the original DGP model is the existence of 
ghost-like excitations [23]. We do not endeavor to discuss such ghost excitation in the present study 
since our aim is the study of the cosmological implications of the model presented in this work. For 
a comprehensive review of the existence of 4D ghosts on the self-accelerating branch of solutions in 
DGP models see [24]. 

In this paper we consider the dynamics of a scalar field existing in the bulk and study the cosmo- 
logical implications that such a scenario would entail in a DGP brane-world. We do this by studying 
the effective field equations on the 3 -f- 1 dimensional DGP brane which is assumed to be rigidly 
located on the orbifold symmetry along the 5th dimension. In section 3 we consider the weak field 
approximation and show that this theory predicts a modification to the mass appearing in the grav- 
itational potential in the form of an additional mass term related to the effective intrinsic curvature 
and that such a mass can be interpreted as dark matter. Finally, in section 4, we discuss the ensuing 
cosmology, showing that the model predicts an accelerated expanding universe without introducing 
dark energy and is consistent with the present observational bounds on the value of the Brans-Dicke 
parameter uj. Conclusions are drawn in the last section. 



2 Brans-Dicke Brane in DGP Scenario 

Let us start by considering a bulk scalar field non-minimally coupled to gravity in a 5-dimensional 
DGP brane world. Such a scalar field in the five dimensional theory can be viewed as a dilaton which is 
purely an outcome of dimensional reduction from some higher dimensional theory to a 5-dimensional 
space-time [25] . The action for our model can be written as 

^5 = ttV / d^^V^ i<t^T^ - '^dA<pd^^(p] + TTI / d^x^R + [ cfx^Cm , (1) 

2«;^5) J \ (f) J 2lJ,^ J brane J 

where TZ is the Ricci scalar associated with the 5-dimensional space-time metric g^g , 4> is a scalar 
field which we shall call the Brans-Dicke (BD) field, w is a dimensionless coupling constant which 
determines the coupling between gravity and the BD scalar field. Similarly, the second term is 
the Einstein-Hilbert action for the induced metric q^i, on the brane where R is the relevant scalar 
curvature, /jp' = 87rG(4) and Cm represents the Lagrangian for the matter fields. Latin indices denote 
5-dimensional components {A,B = 0, . . . , 5) and for convenience we choose k^^^ = 87rG(5) = 1. The 
induced metric g^,/ is defined as usual from the bulk metric g^g by 

Q^.u = SiS^g^g. (2) 
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The variational derivative of the action equation (1) with respect to yields the field equations 

1 



Gab = T^AB — 7,9 ab'^ = 



1 

2' 



T% + Tf£- + TAB 



where 



and 



-AB 



c 



9ae 



Ttcurv 
AB 



R 



(3) 



(4) 



(5) 



Note that T^£'" is the contribution coming from the scalar curvature of the brane. The equation of 
motion for the scalar field ^ is given by 



3a; + 4 



(6) 



where T = c is the trace of the energy-momentum tensor of the matter content of the 5-dimensional 
space-time. Notice the factor 'iu+A in the denominator on the right hand side of the BD field equation 
instead of the familiar 20; + 3 in the 4-dimensional case [26]. This is determined by requiring the 
validity of the equivalence principle in our setup, see [27] for a discussion of this topic in the context 
of 4-dimensional Brans-Dicke theory. 



3 Brans-Dicke DGP brane model in the weak field approximation 
regime 

Although equations (3) and (6) represent the more usual or standard form of Brans-Dicke equations 
in the DGP model but to consider the weak field approximation of the model we shall work in the 
Einstein frame [28, 29] given by 



Gab = Q 



1 

'■,B - ^9ab^ 



:,C 



I rj^curv 
-r J-AB 



(7) 



nhi{g4>) 



■ [f + r^n> 



3a; + 4^ 

which is obtained from (3) and (6) by making the transformation 

9aB =Q<t>9AB^ 



(8) 



(9) 



Tab = Q ^<t> ^Tab 



and 



T~'curv — 1 jl — IrTicuru 

AB — y 9 4 



^AB ) 



(10) 



where Q is an arbitrary constant and the tilde on GaB: Tcu™ j, j^gg^^g ^^^^^ these quantities 
arc calculated using the conformal metric g^^g . 

In the weak field approximation of Brans-Dicke theory, in addition to 



we must also assume that 



9aB ~ VaB "I" ^AB ' 



(11) 

(12) 
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where (p^^^ = (p^^\x'^, y) is a first-order term in the energy density and ^Ll ^ 1. Taking into account 



(12) and setting Q = the transformation equations (9) and (10) become 

9aB ~ VaB + ^AB 1 



(13) 



Tab = (1 - <I>^^^G)Tab = Tab (14) 

and 

AB = (1 - r 'y)TAB = Tab > (15) 

where 

hAB=hAB+^^'^OVAB, (16) 

and only the first-order terms in the mass and curvature densities have been retained. Note that 
j^curv |-^gg ^Yie brane. Now, substituting (12) in the field equations (7) and taking equations (13) 
and (14) into account, we get 

GAB = g[TAB + TTr]- (17) 
On the other hand, the scalar field equation (6) becomes 

^ 3w + 4 ^ ^ 

It turns out then that equations (17) are formally identical to the field equations of General Relativity 
with g replacing the Newtonian gravitational constant G(5) , dropping the coefficient Sir. Therefore, 
if g^^g (G(5) ,x'^,y) is a known solution of the Einstein equations in the weak field approximation for a 
given Tab, then the Brans-Dicke solution corresponding to the same Tab will be given in the weak 
field approximation just by taking the inverse of equation (9), that is 

gAB{x'',y) = o~'r'~9AB{o,^'',y) = [i-4>^'Hx^^,y)g] 5^,(a,x^2/), (i9) 

or, equivalently, 

h^g{x^^,y) = h^g{g,x'',y) - <j>('\x^,y)gv^g. (20) 

We therefore conclude that the general problem of finding solutions of Brans-Dicke DGP equations 

of gravity in the weak field approximation may be reduced to solving Einstein field equations for the 
same matter distribution. It should be noted that the Einstein tensor Gab which appears on the left 
hand side of (17) must be calculated in the weak field approximation, i.e., taking as given by 
(13). 

Now, it is well known from equation (17) that the gravitational potential in the weak field ap- 
proximation of DGP expressed in Gaussian normal coordinates in the Einstein frame is given by 
[20] 



U{r) - 



2 2 
cos(2/x^r) cos(2//^r)si(2//^r) -| — sin(2//^r)ci(2/x^r) 

TT TT 



Girr 

where the sine and cosine integrals are defined by the following relations 



(21) 



si(x)= TdC^ and ci(a;) = - (22) 
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Note that this equation is the gravitational potential for the mass density p{f) = MS{f). We must 
now solve equation (18), leading to (p^^^ given by 



(1) 



47r(3a; + 4)r' 

Using equations (16), (21) and (23) gravitational potential in Jordan frame given by 

g(M + M'=«™) 



U{r) = U{r) 



47r(3a; + 4)r ' 



(23) 



(24) 



where it seems as if M*^"™ plays the role of dark matter which has a contribution to the gravitational 
potential. It therefore seems plausible that p'^'^"'^ [f) = M'^'^'^5{f) could be considered as a candidate 
for dark matter. Prom equation (24) we see that the Brans-Dicke DGP model predicts a transition 
scale relating the 4 and 5-dimensional behavior of the gravitational potential in the Jordan frame 



r < 



: U{f) 



4M^M(3a>+4)+6g(M+M'^"'-") _ ti'^M 
247r(3a)+4)r- 6-ivr 



+ 



\''DGP J \'-DGP J 



DGP 



■ U{r) = - 



g(M+M 



Here ^^.^p = ^ and 7 ~ 0.577 is Euler's constant 



4 Cosmological considerations 

Before wc discuss the energy-momentum tensor, let us define the five dimensional metric which has 
the following form 

ds^ = g^^dx^dx^ = q^,u{x^, y) dx^'dx'' + b^ix^", y) dy^ , (25) 

where /it, = 0, . . . , 3 and y is the coordinate associated with the fifth dimension and we will adopt 
a brane-based approach where the brane is the hypersurface defined by y = 0. We also assume an 
orbifold symmetry along the fifth direction y —y. As we shall see in the coming sections, this will 
help us to simplify our calculations. Next we define the energy-momentum tensor 

T^B = r^slbulk + T'^Blbrane ; (26) 

where the subscripts "brane" and "bulk" refer to the corresponding energy-momentum tensors. For 

simplicity we assume that the bulk is devoid of matter other than the Brans-Dicke scalar field. The 
brane matter field is held at y = with the following energy momentum tensor 

r^B|brane = ^diag(-p,p,p,p,0) and rAB|buik = diag(0,0,0,0,0). (27) 

The above expressions are written assuming that the brane is thin and that the bulk is empty. 

Since we are interested in exploring the spatially flat cosmology {k = 0), we consider a 5- 
dimensional flat metric ansatz of the following form 

ds'^ = -r?{T, y)dT^ + a^ir, y)5ijdx'dx^ + b^ir, y)dy'^ , (28) 

where i,j = 1,2,3. With this metric we are now able to write the equations of motion. The (0,0) 
component reads 

a \a I o^\a a \a o J J 



00 



I rpcurv 
"T J-f 



'-0Q 



(29) 
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where 



Too = I 3- + 7 - -T 



a b 2 



(30) 



and 



The components are given by 



rpcurv 

-'oo 



(Cl 



(31) 




n 



n 



where 



n 



2— + — + - - + 2- -- - + 2- 



n 



n 



n 



a h\ b 
- + 2- +- 
a n I 



-2^ + ^ 
a n 



-f-V 



Tf -L T J- T*^" 
ij ' ' ij 



(32) 




and 



d 6 n u; (/) 

2- H h -- 

a b n 2 (j) 



rpcurv 



The (0, 5) component takes the form 

„ I an' ba' 

3 + 7- 

Van b a 



+ 



a' b' n' oj (j)' 

2 1 h - — 

a b n 2 6 



ah d 

+ 2- 

an a 



■■05 ' 



(33) 



(34) 



(35) 



where 



Finally, for the (5, 5) component one has 









(a 






a 


\a 



+ 



n 



where 



^d h uj (j) 
a n 2 (f) 



a, n 
a n 



The equation of motion for the Brans-Dicke field reads 



n\2 r 
b) 



55 



„ a' n' CO S' 
3 I i_ 

a n 2 (f) 



(36) 



(37) 



(38) 



a b h 
3- + --- 
a b n 



n\2 
b) 



a b n 



= —n 



2 (T + T™*^^ 
3a; + 4 



(39) 



where a dot represents the time derivative with respect to r and the prime corresponds to derivatives 
with respect to y. Note that in the above derivation we have assumed (f) = (f){T,y). We make 
the assumption that the metric and the BD field are continuous across the brane localized at y = 0. 
However, their derivatives can be discontinuous at the brane position in the y direction. This suggests 
the second derivatives of the scale factor and the BD field will have a Dirac delta function associated 
with the position of the brane. Since the matter is localized on the brane it will introduce a delta 



6 



function in the Einstein equations which will be matched by the distributional part of the second 
derivatives of the scale factor and the BD field. 

Using the Einstein equations it is possible to find out the jump conditions for a and n by matching 
the Dirac delta functions appearing on the left-hand side of the Einstein equations to the ones coming 
from the energy-momentum tensor equation (26). For the BD field one has to use equation (39) to 
evaluate the jump conditions. We therefore find 



n 



no 60 



(3a; + 
1 



{3u; + 4)4>, 
2 



p + p™"' + (a; + l)(/9 + /9™"')J , 
3(w + l)(p + p™'^) + (2w + 3)(p + 



(40) 
(41) 

(42) 



where 



P 



IJ?nl 



-2- - 2 

2 



tin Oin 



(43) 
(44) 



7 = l(3u;^-l) and 7c = ^(3u;c - 1), (45) 

with Wc = pcurv and the subscript stands for the brane at y = 0. The first two conditions, equations 
(40) and (41), are equivalent to Israel's junction conditions in general relativity [30] (see [6] for a 
discussion of its application in the context of brane world). It is important to note that the above 
jump conditions at y = depend on the energy density and the pressure component of the brane 
world and induced curvature on the brane. Interestingly, for the radiation dominated phase on the 
brane, p = 3p, the jump condition for cj) does not vanish and is proportional to the energy density 
and pressure component of the induced curvature on the brane. 

Taking the jump of the (0, 5) component of the Einstein equation and substituting equations (40 
and 41) we get the continuity equation for the matter on the brane 

/9 + 3(p + p)^ = 0. (46) 

Equation (46) shows that the energy content of the brane is still conserved in this scenario which 
seems to be at odds with recent results obtained independently by several authors who conclude that 
the presence of a dilaton field in the bulk will lead to a non-trivial coupling with the matter on the 
brane which, from the point of view of an observer living on the brane, would be seen as matter 
leaking from the brane [12, 31, 32, 33, 34]. Our situation is different because the coupling between 
the BD field and ordinary matter on the brane was chosen so as to satisfy the equivalence principle, 
as was mentioned above. Had we chosen the coupling in (39) to have the usual 4-dimensional value 
(2a;-|-3)~^, we would have ended up with a situation where the conservation equation (46) would not 
hold and energy could leak from the brane. 

Taking the mean value of the (5, 5) component of Einstein's equations we can now obtain a 
Friedmann type equation on the brane by following a very similar procedure to the one introduced 
in [6]. Using the fact that due to the orbifold symmetry y —y, the mean value of any of the 
quantities a, n 01 (f) should be zero, we can discard all the terms involving mean values in the average 
of the (5, 5) component of the Einstein equations. The equation so obtained will still involve a term 
containing (p, but we can use the mean value of the BD field equation (39) to write this in terms of a, 
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b and n, and their derivatives. After a somewhat lengthy calculation we obtain the first Priedmann 
type equation on the brane 




1 a' 



6a; 



2a4^2 L(3^ + 4) 

12. 



+ 



+ 1 



270^ 



„• 2 



1 



^- 2 



M2(3w + 4)a2(/)2 2^2(30; + 4) a2 02 
27a; a„ 1 



UJ 



+ 



+ 



/x2(3a; + 4)2ao02 2/x2(3a; + 4)2 ao(/)2 



(47) 



2/x4(3a; + 4) a3<^2 L(3a; + 4) 



-2 -3a; 



+ 



2/x4(3a; + 4)a2</)2 [(3a; + 4) 



6u 



+ 1 



4(3a; + 4)2<^2 



3p^ + 2pp+^ 



2 2 



4(3a; + 4)(^2 

Upon using the components (5, 5) and (0, 5) we obtain the other Priedmann type equation 



(3a; + 4) 



2^2 




1 



dT 



(48) 



dT 

3a4(3a; + 4) ' ^ 



p*"* + (a; + 1)^*"* a^T't'^o 



y=0 



y=0 



While deriving the above equations we have assumed that from the point of view of the brane 
observer, the extra dimension is static, that is 6 = 60. We have also fixed the time in such a way that 
no = 1, corresponding to the usual choice of time in conventional cosmology. It is interesting to note 
that by taking the limit a; — 00 on the right-hand side of equation (47), we obtain exactly the same 
expression as in the DGP model [20]. Note that in this model a; is a coupling constant, but after 
inducing equations on the brane it can be interpreted as the induced BD parameter on the brane. 

Taking the mean value of the BD field equation we obtain an equation of motion for on the 
brane 



^ + 3 



+ 



3a; 



Oo^o /^^(3w + 4)202 
3a; 



6 at, 3 d„ 



(49) 



/x2(3a; + 4)202 «2 



3 a' ^ 

27 + —^ 



a;7 



(3a; + 4)202 • 



Note that in order to obtain equation (49) we also have to assume that the non-distributional part 
of 0" vanishes, otherwise, a term involving 0" will appear in the BD field equation. As we shall 
see in the next section it is possible to obtain cosmologically interesting solutions which verify this 
condition. Using equations (47), (48) and (49) we can determine the cosmology of the brane at y = 0. 



5 Cosmological solutions 

It is well know that the original DGP brane cosmology model leads to self accelerating solutions [17]. 
In this section, we start by consider the self accelerating solutions in our model in the vacuum sector 
i.e p = p = 0. Let us take the brane metric, q^i,, as spatially flat in the 4D de Sitter geometry and 
write 

qi^^dx^'dx" = -dr"^ + e^^^dx^. (50) 
Thus, in order to obtain the self accelerating solutions, we make the ansatz 

(P^{t) = aa^ir), (51) 
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where a is an arbitrary positive constant. Using equations (48), (50) and (51) one finds 



H = ei,\3aj + A)^,^'^^, (52) 

where e = ±1. Note that for e = +1, or the self-accelerating branch, we will have the self accelerating 
solutions in our model. 

We now obtain solutions for a and <p which will allow us to discuss the cosmological implications 
in our setup. Let us then assume that both the total energy density p and pressure p on the brane 
consist of two parts 

p = X + g and p = —A + p, (53) 

where A, g and p arc the tension, the usual cosmological energy density and pressure in the matter 
frame respectively. In what follows we concentrate on the case ^ = p = 0, i.e. the vacuum solution. 
However, by retaining a non-zero effective tension on the brane we are actually taking the brane 
effects into account. For simplicity we take the form of the tension, A, as follows [35] 

A = Xc(l>'', (54) 

where Ac is a constant. As before, we assume that the extra space-like dimension is stabilized (b is 
constant) and choose time such that uq = 1. We now look for a power law solution for the scale 
factor. Substituting the ansdtze 

aQ{T)ocT'^ and 4'oi''') = ~2' 

into equations (47) and (49) one obtains 

m = —1, (56) 

and from equation (49), n is obtained as the solution of the following algebraic equation in terms of 
Ac, fJp and uj 



36a; 4 36a; o n 
rn — — —T^n + n 



(3a; + 4)2'" (3a; + 4)2'" [(3a; + 4)2 p^{3iv + 4)\ 



9a; 24Aca; 



(57) 



— n 



12Aca; 

3 + 



M2(3a; + 4)2 



4A2a; 

+ M^(3a; + 4)2 + ^ " 



This algebraic equation has four explicit real solutions for n in terms of Ac, /i^ and lo. 

Equation (55) shows that there is a possibility of having an accelerated expanding universe for 
some choices of a;, /x^ and Ac. The deceleration parameter on the brane as a function of a;, p^ and Ac 
is given by 

g(a;,Ac,/x2) = -^ = -^, (58) 

where n is one of the four roots of equation (57). A glance at equation (58) reveals that using 
the condition for acceleration, that is q{u>, Ac, p'^) < 0, leads to n > 1 from which, using definition 

= — 1 — 3^ for the effective quintessence, we find w^s < —1/3. Figure 1 shows the behavior 

of the deceleration parameter, q, as a function of lo for % ~ 1 and n, taken as a root of equation 
(57). Therefore, each graph in this figure corresponds to one of the roots of the fourth order algebraic 
equation (57). As can be seen, the graphs in the first row in Figure 1 show that for positive and 
negative values of cu, we have an accelerating universe. The graph on the top left hand corner is 
particularly interesting since it shows that for oj — > ±00 we have g — ^ — 1, that is the universe finally 
approaches the eternal de Sitter phase. It has been claimed that the value of |a;| in 4D should be 
large (> 40000) if Brans-Dicke theory is to be consistent with the astronomical observations [36]. 
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Figure 1: The behavior of q{uj) as a function of w for ^ ~ 1. Each figure represents the deceleration parameter for 
one of the four roots of equation (57). 

Therefore, the graph mentioned above seems to be in good agreement with the observational data. 
The graphs in the second row in figure 1 predict a decelerating universe which is not consistent with 
the present observational data. However, the model can be used to interpolate back the deceleration 
parameter to an earlier epoch to yield a decelerating universe which is required in order to explain 
processes like nucleosynthesis. 

6 Conclusions 

We have derived the effective four dimensional field equations in a DGP brane-world setting where 
a BD field is assumed to be present in the bulk. We considered the weak field approximation in our 
model which resulted in the modification of the mass appearing in the gravitational potential by the 
addition of the mass of the effective intrinsic curvature that appears on the brane. In our scenario, the 
conservation equation for the matter fields confined to the brane still holds in spite of the existence 
of a BD field in the bulk. This is due to the fact that our calculations were done in the Jordan frame. 
Finally we obtained the modified Friedmann equation on the brane and showed that an accelerated 
expanding universe results for certain choices of the parameters which is in good agreement with the 
present bounds on the value of ui in Brans-Dicke theory. 
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